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The time evolution of random variables with Levy statistics has the ability to develop jumps, 
displaying very different behaviors from continuously fluctuating cases. Such patterns appear in 
an ever broadening range of examples including random lasers, non-Gaussian kinetics or foraging 
strategies. The penalizing or reinforcing effect of the environment, however, has been little explored 
so far. We report a new phenomenon which manifests as a qualitative transition in the spatial 
decay behavior of the stationary measure of a jump process under an external potential, occurring 
on a combined change in the characteristics of the process and the lowest eigenvalue resulting from 
the effect of the potential. This also provides insight into the fundamental question of what is the 
mechanism of the spatial decay of a ground state. 
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I. LEVY STATISTICS IN PHYSICS 

Jump processes and related objects (Levy flights, con¬ 
tinuous time random walks, heavy-tailed probability 
distributions, non-local operators, fractional differential 
equations) are increasingly studied in mathematics and 
the natural sciences due to their interest in modelling a 
multitude of phenomena, ranging from physics, through 
chemistry and biolo^, to social science. Surveys with a 
wide scope include [ll-@. 

For a long time, Brownian motion played a central 
role in these modelling efforts, which has complex but 
continuous sample paths, and other appealing statistical 
features such as moments of any order and a scaling prop¬ 
erty. These ideas are now significantly extended by jump 
Levy and Levy-type (e.g., Feller) processes, whose appli¬ 
cation allows to take account of further features such as 
sudden changes, spiky time series, heavy-tailed distribu¬ 
tions, and are increasingly recognized to provide life-like 
additions and refinements to the previously used concepts 
and tools, leading to a new paradigm in scientific mod¬ 
elling. Due to the complexity of the problems, and the 
difficulties and subtleties accompanying them (possible 
non-existence of statistical parameters such as any mo¬ 
ments, non-locality, long range memory effects) there is a 
strong need to have an ongoing communication between 
mathematics and physics. To address these challenges, 
a concurrent mathematical research on the theoretical 
foundations of these new techniques reached a high level 
of development [Mi- 

Recall that by a Levy process one understands a 
stochastic process {Xt)t>o with independent increments 
in the sense that if0<r<s<t are any distinct time 
points, then the random displacement Xt—X^ is stochas¬ 


tically independent from Xg — Xr- (Here Xt is to be 
understood as a position in physical or a value in phase 
space, taken at time t at random with some distribu¬ 
tion.) Furthermore, the increments are also stationary 
in the sense that the probability distribution of the ran¬ 
dom displacement Xt — Xg only depends on the difference 
t — s, which means that displacements over equally long 
time intervals are identically distributed. Thirdly, with 
probability one the paths t of a Levy process are 

continuous from the right with left limits, i.e., Xt+ = Xt 
but it is not necessarily the case that Xt- = Xt, so a 
jump at t may occur. (Here t— means the moment of 
time infinitesimally preceding t, and t-\- means the mo¬ 
ment infinitesimally soon after t). In fact, as it turns out, 
Levy processes having jump discontinuities are the rule 
rather than the exception, and Brownian motion is rather 
the exceptional Levy process which has continuous paths 
with probability one. Due to this special role of Brownian 
motion, the paths of a Levy process (Xt)t>o can in gen¬ 
eral be decomposed into four independent components so 
that Xt = X^^ + Xf ^ -b Xf ^ -b Xt^'‘\ where X^^ is the 
drift component of the form bt (with drift coefficient b) 
having the effect of building in a non-random tendency 
into the full process, Xj ^ is the Brownian component, 
Xt is the small-jump component (with jump size up 
to 1), and Xj^^ is the large-jump component (with jump 
size larger than 1). This decomposition implies that the 
paths of a general Levy process contain stretches of con¬ 
tinuous but very rugged curves, interspersed by jumps of 
any size and direction occurring at random times. 

The frequency of jumps is described by the so-called 
Levy measure v{z)dz, constructed in the following way. 
A jump at time t occurs when AXj = Xt — Xt- is non- 
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zero, thus the total number of jumps up to time t falling 
in a d-dimensional box E can be counted by N(t,E) = 
#{s £ t] : AXs £ E}. By general theory, it is then 
known [sj that the random number N{t, E) follows a 
Poisson distribution whose mean is tv{E) = t v{z)dz, 
thus v{z), z € E, measures the frequency of jumps of size 
in E. Moreover, the Levy measure satisfies the integra- 
bility condition min{l, |zp}z/(z)(iz < oo, which 

means that the total activity v{z)dz of the large 

jumps is finite. Note that this does not rule out the 
possibility that the small-jump activity /o<;|j,|<i v{z)dz 
is infinite. 

All of the independent components of the full process 
contribute by a separate term into the exponent 'L of 
the Fourier transform (characteristic function) of AT^, i.e., 
jgo 1 ‘giy.Xtj _ which is given by the well-known 

Levy-Khintchine formula 

«'(y) = -ib-y+^y-Ay (1) 

+ [ + iy ■ zln^i<i}{z))v{z)dz. 


Here E'^j...] = J ...dF^ denotes expectation (average) 
with respect to the paths of the process starting in 
a; £ distributed by the path measure and 1e{z) 
is the indicator function of set E, i.e., it equals 1 if z £ i? 
and 0 ii z ^ E. The vector b is the drift coefficient, the 
matrix A is the diffusion matrix, and i'{z)dz is the Levy 
measure. 

Levy processes are widely encountered in physics. A 
first observation of the utility of Brownian motion de¬ 
rives from Feynman’s path integration approach to quan¬ 
tum mechanics. Making a Wick rotation t i— >■ it, the 
Schrodinger equation f9t^(a;, t) = {—L+V)'ip{x,t), where 
L = iA is half the Laplacian on L^(R‘’*) and a system 
of units is adopted in which h = 1, is transformed into 
a diffusion equation with a dissipation function given by 
V. Then, with initial wave-functio n ibix, 0 ) = (j^ix), the 
solution can be represented as (see [lj-[ig for details) 


tp{x,t) = (e ^+'^'i(j)){x) = [ T{x,0-,y,t)cj){y)dy 


= E’" 


e-/o 


( 2 ) 


where E’^ denotes averaging with respect to the paths of 
Brownian motion {Bt)t>o- The right hand side thus gives 
a stochastic representation of the time-evolution semi¬ 
group Tt = with propagator T{x,0;y,t), and 

the problem of studying the solutions of the Schrbdinger 
equation is now a problem of statistical mechanics in 
which a random mover is performing Brownian motion 
in the potential landscape V. 

A similar Feynman-Kac type representation holds for 
the semi-relativistic Hamiltonian with kinetic term —L — 
(—A + — TO, involving the fractional Laplacian of 

order ^ and where to > 0 is the rest mass of the par¬ 
ticle |17l - lli^ . Then the corresponding semi-relativistic 


Schrbdinger equation transforms into an anomalous (i.e., 
non-Gaussian) diffusion equation, and Brownian mo¬ 
tion is replaced by a jump Levy process, specifically, 
a Cauchy process with or without mass. When the 
mass is non-zero, the Levy measure behaves like e{z) ^ 
e-"*bl|2;|-^-i(l-I- \z\Y^'^, i.e., the frequency of jumps de¬ 
creases in leading order exponentially with their size. 
When the mass is zero, v{z) ^ i.e., the jump 

size distribution has polynomially heavy tails. (Here and 
below ^ denotes large \z\ asymptotic comparability). 


Quantities obeying Levy statistics occur also in many 
other models or experiments in physics. The Cauchy 
process is the time evolution of stable random variables 
with index 0 = 1 , while other values of a have further 
relevance. The Breit-Wigner formula for high-energy res¬ 
onances is another example of occurrence of a symmetric 
Cauchy distribution Particles travelling between a 
point source and a detection plane produce an asymmet¬ 
ric Cauchy distribution, uniformly distributed stars in 
space give rise to a gravitational field with stable dis¬ 
tribution of index a = 1.5, the size of large polymer¬ 
ized molecules a distribution with the same index, the 
flight time of particles trapped in the vortices of a flow 
is described by a = 1.3, and so on [2lj. In laser cool¬ 
ing technologies, based on making atoms accrue and get 
trapped in low-momentum regions of the phase space, the 
trapping time is a stable random variable with a = 1.5, 
and the recycle time (the time needed to return to the 
trapping reg ion), is a quantity with a = 1.25 [l^, [2^, 
see also [^. Laser optics also provides Levy statistical 
behaviors involving exponentially light-tailed jump mea¬ 
sures. In it has been reported that a medium of Levy 
scatterers has been experimentally manufactured to pro¬ 
duce a “Levy glass”, i.e., an optical material in which 
light travels by following the rules of a Levy flight. In or¬ 
der to implement finite size requirements, so-called trun¬ 
cated Levy processes have been proposed , which 

set hard cutoffs on jump sizes. A variant of these models 
works with soft cutoffs by exponentially suppressing large 
jumps, involving Levy measures of the form 
[281. 129 I]. In a recent development it has been found that, 
since excitation pulses carry finite energies, in random 
lasers the output intensities can be statistically described 
in terms of such truncated and exponentially tempered 
Levy variables, where parameter fittings yielded the val¬ 
ues (5 = 0.53, P = 0.0054 for unscattered light, and 
other truncations for the scattered light component [30l| . 
Exponentially truncated situations have been observed 
also in sub-diffusion models [^ . Purely exponential and 
Weibull distributions are supplied from extreme value 
phenomena, used in the description and prediction of 
the behavior of e.g. maxima of a time-series, related 
to natural catastrophes, financial crashes, technological 
accidents etc [3^. The examples are further multiplied 
to a great extent by models of anomalous transport the- 
ory [il, i, [Hi! statistical physics [lia chemistry 
[41I . [4^ , biology [43I - I45II , or social systems . 
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II. LEVY MOTION IN AN ENERGY 
LANDSCAPE 

In order to unify the wide range of natural phenomena 
leading to Levy statistics in a single framework, we now 
consider a model of a particle performing a spherically 
symmetric R'^-valued Levy motion {Xt)t>o, with diffu¬ 
sion matrix A = aid (where Id is the identity matrix in 
K^) and Levy intensity of the generic form 

= (3) 

with parameters c > 0, /3, <5 > 0, under the effect of an 
external potential V. Our goal in this paper is to study 
the stationary behavior of such a random mover in func¬ 
tion of the jump measure and the potential. Specifically, 
we focus on the spatial decay properties of the stationary 
measures, and report a qualitative transition in the decay 
rates as the jump measure is varied from having a heavy 
to a light tail. Although in our argument we motivated 
the statistical mechanical problem by the Feynman-Kac 
formula-based approach to non-relativistic and relativis¬ 
tic quantum mechanics, it should be clear from the above 
that there are many anomalous kinetic models using non¬ 
local operators and Levy processes to which our frame¬ 
work can be applied. Currently there are many efforts 
trying to find solutions of non-local equations, which is 
a difficult problem. Precise information on the asymp¬ 
totic behavior of eigenfunctions of non-local operators is 
important in verifying conjectures and the validity of nu¬ 
merical approximations. 

We distinguish the following three categories of jump 
measures: 

1. sub-exponential Levy measures: /3 € [0,1), where 
/3 = 0 corresponds to purely polynomial decay, and 
P S (0,1) corresponds to stretched exponential de¬ 
cay in leading order 

2 . exponential Levy measures: P = 1 

3. super-exponential Levy measures: P > 1. 

These jump types cover a large number of Levy pro¬ 
cesses of interest. Apart from the examples discussed 
above, there are many other processes which fit in one 
of these classes (jump-diffusions, geometric Levy pro¬ 
cesses, a variety of other subordinate Brownian motions, 
Lamperti-type processes, hyperbolic processes etc). Due 
to a widespread occurrence of scale-invariant phenom¬ 
ena and self-similar structures (fractals, ubiquitous power 
law-distributed data from natural and social empirical 
observations etc), and of the theory of regular variation 
in mathematics, there is a rich understanding of the sub¬ 
exponential class to date, while exponential and super¬ 
exponential cases were less explored so far [i^, . An 

interesting coinparative study of these three types has 
been made in [sol . [sijj , analyzing the transition times in 
the Arrhenius-Eyring-Kramers laws of chemical reaction 
rates, revealing the differing escape mechanisms through 


a potential barrier in the weak noise limit, according to 
the specific choices of the driving process. We discuss 
here differing qualitative behaviors occurring around a 
sharp transition point of similar jump processes, at arbi¬ 
trary noise strengths. 

With these entries the Levy-Khintchine formula m re¬ 
duces to 


^(y) = + / {l-cos{y z)y{z)dz, (4) 

^ Jr'^\{o} 


and in the path decomposition described above we have 
a Brownian component with diffusion coefficient a > 0, 
and (small and large) jump components with Levy mea¬ 
sure i'(z)dz, while the assumed symmetry prevents a drift 
term to occur. Spherical symmetry implies z^(z) = i^dz]) 
and thus the process can make a jump of size |z| in any 
direction with equal likelihood. 

Next we consider such a random mover in an energy 
landscape described by a potential V, acting as a mecha¬ 
nism reinforcing or penalizing the random mover to go in 
specific regions of space. We assume that V = E+ — V~ 
is Kato-regular in the sense that its negative part satisfies 


lim sup 



|r-(x,)|ds 


= 0 , 


and the restriction to every ball in of its positive part 
V~^ satisfies the same condition. This condition implies 
by the Markov property of the process that — f* V (Xs)ds 
is exponentially integrable for all t > 0, and thus the 
Feynman-Kac formula holds. On the other hand, it is 
general enough to accommodate potentials which may 
have local singularities and any number of local minima 
or maxima. The crucial fact is that the long time behav¬ 
ior of the process depends on the asymptotic behavior of 
the potential at infinity. We consider the following two 
large classes: 

1 . confining potentials: with lim| 3 ,|_joo V{x) = oo 

2. decaying potentials: with lim| 2 ;|_,,oo V(x) = 0. 

Examples of confining potentials include the quadratic 
(harmonic) and all the even anharmonic potentials, and 
some decaying potentials are potential wells and Yukawa 
potentials (there is a large supply of further examples in 
either class). 

In the model we consider, the potential conditions 
the Levy motion {Xt)t>o turning it into a new pro¬ 
cess (Xt)(>o. Under V correlations are introduced be¬ 
tween different values Xs, Xt, hence this is not a Levy 
process. Notwithstanding this crucial difference, it is 
still a Markov process. By multiplication through the 
Feynman-Kac factor e“-^o the path measure in 

the expectation at the right hand side of ([2]) loses its 
normalization, so Tt is not a stochastic semigroup. To 
determine the stochastic semigroup describing {Xt)t>o, 
consider the ground state (po of the Hamiltonian —L + U, 
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i.e., the unique eigenfunction lying at the lowest eigen¬ 
value Aq of the Hamiltonian, whenever this exists. Note 
that the ground state has no nodes and <^o{x) > 0 for all 
cc S A combination of the eigenvalue equation and 
and division at both sides by ^q{x) gives 

1 = / T{x,0;y,t) ‘^°^^\ dy = f f{x,0;y,t)dP{y), 

jRd f0[X) jRd 


with 


T{x,0;y,t) 


e^°*^T{x,0-y,t) 

<fo{x)(po{y) 


and P given by eq. (I5|) below. Thus the semigroup Tt 
with propagator T{x,0;y,t) > 0 for all x,y S R*^ and 
t > 0, satisfies the normalization Tjl = 1 at all time, and 
so it is a stochastic semigroup with stationary measure 

dP{x) = ip^{x)dx. (5) 

A further computation gives the transition prob¬ 

ability density of the random process (At)t>o from posi¬ 
tion x at time s to a position y at time t to be 


p(x,s| j/,<) 


gAo(t-s) 

y}o{x)ipo{y) 




( 6 ) 


where p(x, s | y, t) is the transition probability density of 
the underlying process {Xt)t>o, and the notation indi¬ 
cates conditional expectation with pin-down conditions 
Xs = X and Xt = y. The path measure of the process 
under V can also be calculated, and we obtain that the 
probability that {Xt)t>o starting from x will be in a re¬ 
gion A in space at time t is 


P“(AteA)= f p{x,0\y,t)‘fl{y)dy. (7) 

JA 

In light of dl]), the right hand side shows that the path 
measure can be interpreted as a Gibbs distribution over 
the paths of the random process, for the energy function 
given by the integral of V. 

The infinitesimal generator of {Xt)t>o is (see also M) 

(Lf){x) =^Af{x)+aVlogipo{x)-Vf{x) 

f ipo{x + z) - (po{x) ^rf ^ / N , 

+ / -- z-V f{x)v{z)dz 

Jo<\z\<l <f 0 (x) 

+ [ {f{x + z)-f{x) 

The expression indicates that the jump measure and drift 
are now position dependent (unlike for a Levy process), 
and the process (At)t>o will have different jump prefer¬ 
ences according to what the value of H at x is. Note that 


the effect of the potential comes in via the ground state 
ipo- Interestingly, for the process constrained by the po¬ 
tential two types of drift appear. The first takes the usual 
form and is related to the diffusive part of the underlying 
process. The second corresponds to the jump part of the 
underlying Levy process and is given by the integral dif¬ 
ference operator. While now the motion is in some sense 
more complex and loses its isotropy, the topology of the 
paths is preserved on applying the potential, in the sense 
that Brownian motion under V keeps having continuous 
paths, and a jump Levy process under V keeps having 
continuous stretches interrupted by random jumps. 


III. STATIONARY BEHAVIOR 

In general, when V = 0 the process {Xt)t>o will dis¬ 
play wild fluctuations, exploring the full space. Brown¬ 
ian motion has no t —>■ oo limiting behavior in the sense 
that liminft_>oo Bt = —oo and limsup(_^oo Bt = oo with 
probability one. However, the law of iterated logarithm 
says that with probability one \Bt\ < \/2t log log t, after 
a sufficiently long time. In other words. Brownian paths 
eventually get localized with a high probability inside a 
curved cone given by the iterated log-profile function. For 
jump Levy processes which have infinite variance (such as 
stable processes), the situation is very different as there 
exists no such profile function and no region in space in 
which paths concentrate on the long run. 

However, when a non-zero potential is applied to 
a Levy process, the long term behavior dramatically 
changes. Paths will now tend to spend long times at the 
local/global minima of the energy landscape since ([6]) 
implies that strays involve an exponentially high price. 
On the long run, the process will tend to settle on a be¬ 
havior described by the stationary measure P given by 
whenever this exists, thus the key to the stationary 
behavior of the process {Xt)t>o is the ground state (po. 
Explicit calculations of ground states for Levy processes 
are known so far for V (x) = x^ and V (x) = x"' . 

The tail behavior of P, which amounts to the asymptotic 
behavior of po as x ^ ±oo, tells then of how strongly 
the random mover is localized in the bulk, in particular, 
around the local and/or global minima of the potential. 

A ground state, and thus a stationary behavior of the 
process conditioned by V, occurs only if the potential 
is suitably chosen. When the potential is confining, a 
ground state always exists. When the potential is decay¬ 
ing, it is in general a hard problem to determine whether 
for a given Levy process and a given potential a ground 
state exists, and it is possible to construct many exam¬ 
ples of jump processes and potentials for which it does 
not [la, [53|- The physical reason why a ground state 
Po forms at all is due to a mechanism allowing the pro¬ 
cess to accumulate a sufficient mean total sojourn time 
in each unit-ball neighborhood of every point in space 
through visiting and revisiting them |58l . l59l | . Clearly, if 
V is confining, then ([5]) implies that far out the paths 
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will be exponentially penalized and a strong centripetal 
effect drives most of the motion to take place near the 
bottom region of the energy landscape, which makes the 
process steady down to a stationary state. If V is decay¬ 
ing, then far out the motion increasingly resembles free 
motion (which has no ground state) and the energetic ef¬ 
fect becomes more delicate so that a ground state forms 
for a sufficiently low-lying ground state eigenvalue only. 

Whenever a ground state does exist, we can derive a 
useful representation by using the eigenvalue equation. 
Since the ground state eigenvalue Ag may be negative 
dependent on the choice of potential, we shift it by 9 so 
that 9 + Xq > 0. Then on multiplying both sides of the 
eigenvalue equation by and integrating with respect 
to time, we get 


ipo{x) = {9 + Ao) 



,-f*{e+v{Xs))d8 


‘foiXt 


dt. 

( 8 ) 


We can manipulate the integral at the right hand side by 
choosing an arbitrary bounded open set D in space, con¬ 
taining the starting point of the process, and considering 
the first exit time 


td = inf{t > 0 : Xt ^ D} 

of the process from D. Using this, the process can be 
stopped at this random time and then allowed to run 
from this time onwards by using its Markov property. 
By splitting up the time integral accordingly, this leads 
to 


ipo{x) = {9 + Ao)x 




= (0 + Ao)E" 




e-/o(®+^(^U)ds 




{rD<oo} 


(poiXt)dt 




It is in our gift how to choose D optimally to study the 
large |a;| behavior of the ground state, and we take it 
to be a unit ball Bi(x) centered in x. Making use of the 
above representation, after some analysis we find that the 
second term provides the main contribution into a lower 
bound of and the first term in an upper bound. 

Putting these estimates together, we obtain 


spent in specific locations of space, the decay of a ground 
state is governed by another mechanism, which deter¬ 
mines how soon the random mover leaves such unit-balls 

far out [^.1^. 

If U is a confining potential with a mild additional 
property, we find through an analysis of the mean exit 
times in m that whenever v is sub-exponential or expo¬ 
nential with 6 > , the ground state behaves like 


^oix) 


v{xy 


This gives a neat account of the separate contributions 
of the potential and the underlying free process, and the 
expression shows that the decay rate of tpQ results from a 
balance between the Levy intensity and the killing effect 
of the potential far out. 

If U is a decaying potential, the situation becomes 
more subtle since now both v and V go to zero and a 
similar balance does not hold. The effect of the potential 
comes in now through the negative ground state eigen¬ 
value Aq it produces (if any). The following behaviors 
occur: 


(a) polynomially decaying Levy intensities: if v{\z\) ^ 
\z\~^ with S > d, then (poix) ^ \x\~^ for all Ag < 0 

(b) stretched-exponentially decaying Levy intensities: if 

i^d^l) ^ with c > 0, /3 S (0,1), ^ > 0, 

then g}o{x) ^ \x\~^ for all Ag < 0 

(c) exponentially decaying Levy intensities: if v{\z\) ^ 
e-chl|0|-'5 with c > 0 and 5 > then there is 
a process-dependent cutoff rjx such that whenever 
|Ao| > T]x, wehave(^o(x) ^ e~'^''^''\x\~^] when there 
is no restriction on the eigenvalue Ag, we have for 
large enough \x\ that for every e > 0 there is a 
Cg > 0 such that 

Mx) > C, max {e-VUol+d^d ; 

moreover, if ^ then ipo{x) > Cv{x) for every 

constant C > 0, at least for large enough |a;| 

(d) super-exponentially decaying Levy intensities: if 

v{\z\) ^ with c > 0, /3 > I and <5 > 0, 

then for every e > 0 there is a > 0 such that 


LP^X) E"' 




-yViX,)ds 


v{x) (9) 


for a large class of Levy processes (Xt)t>g, where tb^(x) 
is the first exit time from this ball by the process start¬ 
ing at X. The factor multiplying v{x) in (jH]) measures the 
mean time spent by {Xt)t>o perturbed by the potential 
U in a radius one neighborhood of x before the first exit 
from this region. This shows that while the creation of a 
ground state depends on the right amount of time being 


(po{x) > 

for large enough |a;| and any Ag < 0 

(e) diffusions: if {Xt)t>o is Brownian motion, V is 
bounded and continuous, then for all 0 < £ < |Ag| 
there is Ce > 0 such that 


^o(^) < 
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From the above behaviors it is seen that the ground 
state decay follows the decay of the Levy intensity for 
sub-exponential cases and Aq which can be arbitrarily 
close to zero. For super-exponential cases the rate of 
decay is slower than that of zz, whatever the magni¬ 
tude of Ao, and this comparison can be further con¬ 
tinued to diffusions to find that (po(x) decays at infin¬ 
ity much slower than p(x, 0 | 0,1), where p(x, sj y,t) = 

^ _ \y-x\^ 

(2Tr(t — s))“ 2 e 2 {t-s) is i;] 2 e Gaussian transition proba¬ 
bility density which can be used instead of v. The ex¬ 
ponential case is an interesting and delicate borderline 
situation, marking a transition line in the decay pat¬ 
terns. In this case also the magnitude of the ground state 
eigenvalue and the sub-leading order polynomial factors 
in ^{x) play a role. Whenever 6 > the absolute value 
of Ao determines what the decay rate of po is: if |Ao| is 
large enough, then the ground state follows the decay of 
the Levy intensity, if |Ao| is too small, then 
may become larger than \x\~^ and thus the decay of 
the ground state gets slower than that of the Levy inten¬ 
sity. As soon as (5 < the decay becomes invariably 
slower than that of v, no matter the value of |Ao|. 

These results are obtained through an analysis of spe¬ 
cific preferred jump scenarios and the processes’ response 
to perturbation by the potential. A key observation is 
that the sub-exponential processes (with S > d ioi poly¬ 
nomial and 5 > 0 for stretched exponential v) satisfy 
the property that events in which two arbitrary points 
are linked through a sequence of any length of successive 
large jumps via intermediary points, are stochastically 
dominated by events consisting of a single direct large 
jump between the given two points. This jump-paring 
property has been explored in [^, |6l|, It fails to 

hold for super-exponential processes, and there is a di¬ 
vide around the critical exponent 5 = for exponen¬ 
tial processes, above which it does hold, and at and below 
which it breaks down. When the jump-paring property 
fails to hold, the ground state necessarily decays slower 
than the Levy intensity. The process-dependent cutoff 
rjx appearing in (c) above relates to how much energy 
needs to be fed to the process to make it “bend” and 
follow the decay of the intensity of jumps. 


IV. TRANSITION IN THE DECAY RATES 

The regime change occurring at exponentially decaying 
Levy intensities can be explained by a switch-over in the 
mechanism which makes the long jumps of the process 
occur. This can be seen by separating the large-jump 
component {Xt)t>o from the full process {Xt)t>o- Its 
Levy intensity can be obtained by restriction of the Levy 
intensity of the full process to large jumps, i.e., Va^{z) = 
z/(z)]{|z|>i}. Using this, the jump-paring condition is 


formulated as 

{Voo * i'og){z - x) = / Vooiy - x)voc,{z - y)dy 

jR'* 

< Cv^{z-x), 

where the star indicates convolution and C > 0 is a con¬ 
stant. This expresses the property that a large jump from 
X io y and then on from y to 2 ; is less preferred than a 
direct large jump from x to z. Continued inductively, it 
follows that the n-fold convolution so 

all multiple large jumps are stochastically dominated by 
single large jumps. 

Since by the general integrability condition on the Levy 
intensity stated above v is integrable with respect to the 
large jumps, Voo can be normalized to the probability 
measure Voo{z)dz = i>oo{z)dz/ ^^^^y^Voo[z)dz, which de¬ 
scribes the distribution of large jumps of size \z\. De¬ 
note the successive jumps of (Ai)t>o by Ji, J 2 , •/a, 
of sizes I Ji I, IJ 2 1,1 Ta I,... They are mutually independent 
random variables described by the probability measure 
P{Ji G dz) = Voo{z)dz. Whenever {Xt)t>o is a Levy 
process with sub-exponentially decaying v, we have for 
large enough r that 

Pi\Ji + J2\>r)^2Pi\Ji\>r). 

A straightforward calculation using this property gives 
P(max{| Ji|, IJ 2 I} > r) Ri P(| Ji -k J 2 I > r) 
as well as 

:p(^<|Ji|<r, ^<|J 2 |<r) = o(P(|Ji|>r)). 

The implication is that a displacement in the random 
mover’s position resulting from a double jump Ji,J 2 
is larger than a large value r when either |Ji| > r or 
IJ 2 I > r, and it is much less likely that r is exceeded 
cumulatively by jumps having comparable size. This 
property extends to arbitrarily long sequences of jumps. 
When (Xt)t>o has an exponentially decaying v, this no 
longer holds and we have (with 5 = 0) 

P (^ < |Ti| < A ^ < \ J2\ < r) Ri P(|Ji| > r), 

i.e., the same order of magnitude, meaning that now a 
large displacement r is achieved typically cumulatively 
rather than through a single large jump. This can be seen 
as an increase of the capacity of the process to fluctuate 
through multiple smaller jumps rather than exceedingly 
large single jumps, which improves as the tail of Poo (and 
thus of v) is chosen lighter. 

The conclusion one can draw from this transition phe¬ 
nomenon is that there are unexpected features in how 
jump processes respond to even small perturbations by a 
potential. To appreciate this on a specific example, con¬ 
sider a d-dimensional potential well V{x) = —ul{|a;|<a} 
of depth V > 0 and diameter 2a, and a random mover 
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performing Levy motion {Xt)t>o of the above jump cat¬ 
egories, in the landscape of this potential. In this case 

/ V{Xs)ds = -v l{\Xs\<a}ds =-vUt{a), 

Jo Jo 

where Ut{a) is the total time the process spends in the 
potential well within the time interval [0,t]. Since in 
d > 3 every Levy process is transient and the mean to¬ 
tal sojourn time E“[17oo(a)] is finite, a ground state will 
form only for sufficiently large depths v > v* oi the po¬ 
tential well, and the threshold value v* depends on the 
process. When d < 2, then for some processes a ground 
state forms for arbitrarily small non-zero v (e.g., in the 
case of Brownian motion), while for others it may de¬ 
pend on further details (e.g., for a spherically symmetric 
a-stable process there is a ground state for every n > 0 if 
d = 1 and 1 < a < 2, but only for large enough v when 
0 < a < 1, while for d = 2 for every a a deep enough po¬ 
tential well is needed). When a ground state does exist, 
the processes can be ranked according to their inherent 
ability to form such a ground state, and we find that 
Brownian motion and processes with super-exponential 
Levy measures form a ground state easiest, and polyno- 
mially heavy-tailed processes the hardest [^ . 

Assuming that a ground state exists, consider now ran¬ 
dom motion according to a stable process, starting from 
the origin, i.e., the middle of the potential well. Suppose 
that after some time the mover gets to a long distance 
r from the origin. Due to the dominance of single large 
jumps, it is likely that it will get there by a single jump of 
comparable size to r, without “wasting” too much time 
with intermediary visits elsewhere. Before the next large 
jump occurs, the process performs smaller fluctuations 
more locally. Since there is an energy gain in returning 
to the region inside the potential well, this will happen 
with a good chance, however, just because another large 
jump is due soon, the mover will not stay around zero for 
long. This back and forth motion repeats then at a rate 
proportional to the jump intensity, and so the process 
develops a stationary distribution tracking the curve of 
v, with their tail behaviors coinciding exactly. A simi¬ 
lar pattern occurs for all motions of sub-exponential and 


jump-paring exponential type (i.e., (5 > ^^). 

Next consider processes with exponentially decaying 
jump intensity for which the jump-paring property no 
longer holds (<5 < ^^). Now the paths from a long 
distance r can cluster around the origin less efficiently 
and once they are back, they spend comparatively large 
amounts of time in its neighborhood, and re-entries at 
shorter gaps build up “backlogs” in the decay-events so 
that m does not hold any longer. Thus their ground 
states will decay possibly much slower than iz. The same 
pattern occurs for super-exponentially light jumps and 
continuously fluctuating processes (diffusions), which 
have strong limitations to make long trips in short times. 
The accumulation of these backlog events makes a sud¬ 
den qualitative change to occur, which can be seen in the 
polynomial order sub-leading term at a critical exponent, 
which we can identify explicitly to be 

The property that the ground state follows the decay 
regime of the Levy intensity is a combined effect of the 
jump-paring property favoring single large jumps over 
multiple large jumps, and the path concentrating effect 
around the bottom of the well due to the ground state 
energy Aq. The efficiency of this mechanism depends on 
how low-lying Aq is. The above results show that when 
the process typically evolves through large jumps, a small 
energetic effect suffices for securing an efficient concen¬ 
tration. As a result, the hierarchy of good stationary 
behavior is then that the ground state of a jump pro¬ 
cess in a deep enough potential well decays like v, decays 
slower than v if the well is not deep enough, and the 
ground state may even cease to exist if the well is too 
shallow. 
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